Laser induced magnetization switching in films with perpendicular anisotropy: a 
comparison between measurements and a multi-macrospin model 
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Thermally-assisted ultra-fast magnetization reversal in a DC magnetic field for magnetic mul- 
tilayer thin films with perpendicular anisotropy has been investigated in the time domain using 
femtosecond laser heating. The experiment is set-up as an optically pumped stroboscopic Time 
Resolved Magneto-Optical Kerr Effect magnetometer, ft is observed that a modest laser fluence 
of about 0.3 mJ/cm^ induces switching of the magnetization in an applied field much less than 
the DC coercivity (0.8 T) on the sub-nanosecond time-scale. This switching was thermally-assisted 
by the energy from the femtosecond pump-pulse. The experimental results are compared with a 
model based on the Landau Lifschitz Bloch equation. The comparison supports a description of the 
reversal process as an ultra-fast demagnetization and partial recovery followed by slower thermally 
activated switching due to the spin system remaining at an elevated temperature after the heating 
pulse. 

PACS numbers: 75.50.Ss 75.40.Mg 75.40.gb 76.60.es 
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I. INTRODUCTION 

There is currently considerable interest in ultra-fast 
laser-induced magnetization processes. Since the demon- 
stration by Beaurepaire et alX that the magnetiza- 
tion can respond on the picosecond timescale to heat 
pulses produced by femtosecond lasers, a number of 
groups^i^ have studied magnetization processes on this 
timescale. Experiments generally use pump-probe pro- 
cesses in which a high energy laser pulse is used to heat 
the sample and a low energy probe pulse (split from the 
main pulse) is used to monitor the magnetic response us- 
ing the Magneto-Optical Kerr Effect (MOKE). Much of 
this work has investigated the dynamics of the destruc- 
tion and recovery of the magnetization, which can occur 
on the sub-picosecond timescale, although the recovery 
can take an order of magnitude longer due to frustration 
effects among large numbers of nucleation sites at which 
the recovery starts locally^. 

The dynamics of reversal during a pulsed laser experi- 
ment in the presence of an applied field has received less 
attention. Hohlfcld et al.^ investigated the magnetiza- 
tion reversal induced by 100 fs laser pulses in a GdFeCo 
Magneto-Optical recording medium with perpendicular 
anisotropy. They observed an ultra-fast demagnetization 
of the film occurring within the first picosecond; followed 
by a slower recovery, taking nearly a nanosecond, in the 
direction of the applied field as the heat drains from the 
media layer. They analyzed their measurements using 
the Bloch equation and concluded that the reversal pro- 
cess was governed by the nucleation and growth of do- 
mains in the applied field. 

Laser assisted magnetization processes have consid- 
erable potential for future ultra-high density recording 
systems. Essentially, the path to higher densities re- 
quires continuous reduction in the grain volume V of 
the storage medium, necessitating an increase in the 



magneto-crystalline anisotropy energy density K in or- 
der to preserve a sufficiently large value of the parameter 
KV/k-eT (~ 60) to ensure thermal stability of written 
information^. However, the large value of K impacts the 
writability of the information, and some scheme is re- 
quired to overcome this problem. A promising solution is 
Hybrid or Heat Assisted Magnetic Recording (HAMR)i 
in which the medium is heated in order to lower the 
anisotropy and thereby allow information to be written to 
the medium. Since this is a relatively new field the mag- 
netization reversal mechanisms are not well understood. 
Although the work of Hohlfeld et al^ has demonstrated 
thermally activated domain processes in magneto-optical 
media, perpendicular recording requires relatively weakly 
coupled granular media in which domain processes are 
not the dominant reversal mechanism. 



This paper presents time-domain measurements of the 
magnetization reversal process induced by an ultra-fast 
laser pulse in a thin film with perpendicular anisotropy. 
The film was especially designed to have a low Curie tem- 
perature in order to simplify the experiments. We com- 
pare the results with a computational model using the 
Landau-Lifshitz-Bloch (LLB) cquationiS., which is ideally 
suited to simulation of magnetization processes up to and 
beyond Tc and has been shown-- ^ '^^ to give an excellent 
description of the physics of pulsed laser processes. It 
is concluded that the magnetization response consists of 
a fast demagnetization followed by a slower response in 
which the magnetization evolves into the field direction 
by a process of thermally activated transitions over the 
local energy barriers. Our LLB-micromagnetic model is 
shown to give a good description of the physics of the 
reversal process on both timescales. 
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II. METHOD 

The experiment is a stroboscopic pump-probe exper- 
iment using a strong laser pulse to initiate a change in 
the magnetic state of the sample and a weak probe pulse 
to observe the resulting magnetization dynamics via the 
MOKE. The sample is mounted in a spin-stand which 
first moves the magnetic film through a reset field of mag- 
nitude ~ IT for a duration of 0.2 ms which ensures that 
the sample is in a remanent state before exposure to each 
laser pulse. Fig.[l]illustrates the spin-stand arrangement, 
with the insert showing the temporal field profile experi- 
enced by the magnetic film. After resetting to saturation 
the sample then moves into a perpendicular applied field 
from an electromagnet (field range of ±0.52 T). We note 
that the field applied is always lower than the static coer- 
civity of the sample as measured at room temperature by 
a vibrating sample magnetometer. At the center of the 
pole piece a small hole allows optical access to the sample 
at which point the sample is exposed to the pump-pulse. 
The laser pulses arrive at a rate of 1 kHz but the rotation 
of the spin-stand (about 7,000 rpm) ensures that freshly 
saturated film arrives between pulses. 



region probed. The probe approaches the sample surface 
at near normal incidence; this polar MOKE geometry 
yields sensitivity to the out-of-plane component of mag- 
netization. The reflected probe beam is directed into 
an optical bridge detector which uses a WoUaston prism 
to split the beam into two orthogonal polarized compo- 
nents which impinge on a two segment photo-detector. 
By rotating the detector to the angle where the output 
of the two polarization channels (A & B) are approxi- 
mately equal, then the difference between them (A-B) is 
sensitive to small changes of the polarization angle of the 
reflected probe which in turn is proportional to ruz- The 
sum (A-l-B) is sensitive to changes in reflectivity AR, 
which is associated with temperature changes and stress 
waves due to processes such as lattice expansions. In or- 
der to improve immunity to laser drift the pump beam 
was chopped and lock-in ampliflers used to detect the 
sum and difference signals. This makes the measurement 
sensitive to the difference between the state of the sam- 
ple without the pump pulse applied and the state induced 
by the pump-pulse. Therefore the measurements are rel- 
ative and it is difficult to assign an absolute scale to the 
magnetization changes. 
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FIG. 1: Illustration of spin-stand arrangement used to ensure 
the materials magnetization is reset between pump pulses. 

The stroboscopic experiment uses a Libra laser system 
(made by Coherent) which can produce a 1 kHz stream 
of 1 mJ, 150 fs pulses of 800 nm radiation. This beam 
is attenuated, the probe beam is split off and frequency 
doubled to 400 nm. The pump beam is routed around 
an optical delay line with 17 fs resolution over a 1 ns 
range and then focused at normal incidence to a spot 
approximately 800 /im in diameter on the disk surface. 
The probe beam is polarized and then brought to a 400 
fj,Ta focal spot, centered on the pump beam with a power 
level of about l/20th of the pump. Because the probe 
is only half the diameter of the pump it must be noted 
that there will be a temperature distribution across the 



III. EXPERIMENTAL RESULTS 

For characterization of the sample the quasi-static hys- 
teresis loop was measured. Of particular importance is 
the coercivity, which on this timescale is ^ 0.85T, and 
the saturation magnetization Mg, which is equal to 0.35 
xlO^ A/m. The coercivity of course is already greater 
than the field applied during our pulsed field experiment. 
However, it is important to note that the dynamic coer- 
civity on the timescale of the pulsed laser experiment is 
even larger. The intrinsic coercivity (Hq) and the ther- 
mal stability factor {KV/k^T) were measured by making 
a series of time dependent coercivity measurements and 
fitting to Sharrock's law-- . The intrinsic coercivity is re- 
lated to the anisotropy field and is expected to give a 
reasonable estimate of the coercivity at the nanosecond 
timescale. The value of Hq was found to be 1.4T; a fac- 
tor of almost 3 greater than the maximum field available 
from the electromagnet. Separate measurements deter- 
mined the anisotropy if to be 4 x 10^ J/m^, from which 
a grain size of lOnm was estimated. 

Fig. [2] shows a set of time-resolved measurements on 
the sample. Fig. [5] (a) shows the reflectivity data, Fig. [2] 
(b) shows the dynamic magnetic response for zero ap- 
plied held and Fig. [2] (c) shows the response in the pres- 
ence of a reversing field of 0.52 T. The laser pulse energy 
is varied up to 1.14 //J per pulse (which corresponds to 
a fluence of approximately 0.23 mJ/cm^). This value 
is the largest than can safely be applied to this sample 
as energy fluences above about 0.5 mJ/cm^ damage the 
sample. The sample reflectivity data shown in Fig. [2] 
(a) is a probe of the electron/lattice temperature in the 
system. It indicates that the same temperature profile 
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applied field. Recalling that the dynamic coercivity es- 
timated from magnetic measurements is around 1.46T, 
this demonstrates a significant heat-assistance during the 
sub-nanosecond reversal process. 

We now consider in detail the processes involved in 
the magnetization dynamics, which involves three char- 
acteristic timescales as illustrated in Fig. [5] The ini- 
tial phase involves a rapid demagnetization of the sam- 
ple lagging the change in reflectivity by only 50 fs. The 
demagnetization takes about 500 fs, independent of the 
applied field. This is consistent with the normal expec- 
tation of a rapid demagnetization as previously demon- 
strated experimentallyii^i^ and theoretically^ii^. The 
sample then appears to partially recover its magnetiza- 
tion in the original direction, on a timescale indepen- 
dent of the applied field. It is interesting to note that 
the rate of recovery from the demagnetization peak is 
related to the amount of demagnetization achieved - a 
pattern consistent with the calculations of Kazantseva et 
al.^ . Throughout the whole process the only apparent 
field dependence occurs in the longer time-scale dynam- 
ics (20ps-lns). Over this timescale, for the higher laser 
powers, a gradual reversal of the magnetization is seen. 



FIG. 2: Series of measurements in an applied reversing field 
of 0.52T for different laser pulse energies, a) shows the reflec- 
tivity data, b) the rriz component of magnetization, and c) 
shows the response in the presence of a reversing field 



is generated each time and that the temperature scales 
linearly with pulse energy. The reflectivity measurement 
has three distinct features. Within the first 5 ps is a large 
peak having a width of 350 fs, which corresponds to the 
large rise in electron temperature caused by the arrival 
of the 150 fs laser pulse. The electron system then estab- 
lishes thermal equilibrium with the lattice which creates 
the second, rather broader, peak 20 ps later. In reality 
the peak temperature the lattice reaches is much lower 
than that for the electrons; however, it seems that in 
this sample the change in reflectivity is more sensitive to 
the lattice temperature than the electron temperature. 
The waves seen superimposed onto the lattice temper- 
ature peak during the 3-35 ps time frame are probably 
stress pulses launched into the film by the laser heating 
of the surface, which reflect off the interface with the 
glass substrate. The time between peaks is 12 ps, which, 
given the film and interlayer thickness of 25 nm, suggests 
a propagation speed of 4,100 m/s; consistent with the 
speed of sound in a typical metal. The sample appears 
to cool rather slowly as the lattice temperature has only 
fallen to half its peak value after 1 ns. This reflects the 
fact that in this film there is no heat sink to absorb the 
heat. Fig. [5] (b) and (c) show the magnetization dynam- 
ics in applied fields of and -0.52 T respectively. Fig. [2] 
(c) clearly demonstrates heat assisted reversal due to the 
pulse. For pulse energies above ~ 0.5/^J the sample is 
seen to cool with the magnetization aligning along the 
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FIG. 3: Series of measurements of the temporal variation of 
rriz in various applied reversing fields for a constant laser pulse 
energy. 

A similar systematic trend is also exhibited by the tem- 
poral variation of niz at constant laser power in various 
applied fields, as shown in Fig. [3l The initial disappear- 
ance and recovery of the magnetization is similar for all 
applied fields, but with increasing positive field the long- 
term trend is clearly toward positive magnetization. 

The form of the dynamics, involving a trend against 
the field direction followed by a slow reversal into the 
field direction is apparently somewhat counter-intuitive. 
However, the form of the magnetization evolution can 
be explained by consideration of the different timescales 
associated with processes at the atomic and macrospin 
length scales. Within each grain or 'macrospin' the dis- 
appearance and recovery of the magnetization will be 
governed by the longitudinal relaxation time, which is 
temperature dependent but typically of a sub-picosecond 
timescale, which is consistent with the experimental data 
of Fig. [21(b) and (c). However, after this process the sys- 
tem remains at an elevated temperature for around Ins, 
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so there is a possibility of thermally activated magnetiza- 
tion reversal. This will have a characteristic timescale de- 
termined by the macrospin energy barrier, which is low- 
ered by the reduction of the anisotropy constant at ele- 
vated temperature, but typically has values much greater 
than Ips. On this basis we would expect to see a fast re- 
duction and recovery of the magnetization due to atomic 
processes on the picosecond timescale with a slower relax- 
ation into the field direction due to thermally activated 
reversal of the macrospins. In order to test this hypothe- 
sis in relation to the experimental results we have devel- 
oped a computational model of the laser heating process 
based on the LLB equation which has been shown3 to 
give a good description of magnetization processes over 
both characteristic timescales. 



IV. DYNAMIC MODEL OF LASER HEATING 
PROCESS 

We have developed a computational model of laser- 
induced magnetization dynamics of a thin film with 
perpendicular anisotropy. Consistent with experiments 
we assume a granular microstructure for the film. For 
simplicity we assume, in these initial calculations, a 
mono-disperse grain size and anisotropy. Inter-granular 
magneto-static interactions are included, but the grains 
are taken as exchange decoupled. As mentioned previ- 
ously the Landau-Lifschitz Gilbert (LLG) equation can- 
not be used for models of laser heating since it does not 
allow longitudinal fluctuations of the magnetization^i^. 
In the following, we use the LLB equation^° for the ther- 
modynamic simulation of the laser-induced magnetiza- 
tion switching. As described in detail in RefJ^ii^, the 
LLB equation has been shown by comparison with atom- 
istic calculations to give a remarkably good description 
of the physics of ultra-fast high temperature dynamics. 
The LLB equation can be written as 
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[m, X [m, X (H^ff + Cl)]] 



(1) 



Note, that besides the usual precession and relaxation 
terms in the LLG equation (see Ref.J^ for more details), 
the LLB equation contains an additional term which con- 
trols the longitudinal relaxation. Here, rrii is a spin po- 
larization which is not assumed to be of constant length 
and even its equilibrium value mc{T) is temperature de- 
pendent. a|| and a±^ are dimensionless longitudinal and 
transverse damping parameters. 

The LLB equation is valid for finite temperatures and 
even above Tc though the damping parameters and effec- 
tive fields are different below and above Tc. For T < Tc 
the damping parameters are 
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and for T > the damping parameters are equal, 
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In these equations A is a microscopic parameter which 
characterizes the coupling of the individual, atomistic 
spins with the heat bath. 

Thermal fluctuations^^ are included as an additional 
noise term Ci{t) with I =±, ||, (CKO) = and 



2kBT 



(4) 



where i, j denotes lattice sites and v, rj the Cartesian com- 
ponents. Here, A'^ is the volume of the grains and A/f is 
the value of the spontaneous magnetization at zero tem- 
perature. 

The effective fields are VL\g — —^j^, with / the 
free energy density. The total local field is given byiS 
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with the anisotropy field 
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which makes the z axis the easy axis of the model, and 
the dipolar field 
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Note, that the susceptibilities x; si's defined by xi — 
dmi/dHi. At lower temperatures the perpendicular sus- 
ceptibility X-L is related to the anisotropy K via x_l = 
M^ml/{2K)i^. One problem for the application of the 
LLB equation is that one has to know the functions for 
the spontaneous equilibrium magnetization rriciT), the 
perpendicular (x±(T)) and parallel {x\\{T)) susceptibili- 
ties. Here, we use the functions for FePt gained from a 
Langevin dynamics simulation of an atomistic spin model 
as described ini^. In the following, we assume that these 
functions reflect the correct temperature behavior. We 
normalize the perpendicular susceptibility such that its 
value at 300K is consistent with experimentally deter- 
mined values for CoPt {K = 3.94- 10^ i/vo?) determined 
by the rotation method^®. The input functions are shown 
in Fig. m 

The LLB equation is solved numerically by using 
Langevin dynamics simulations as described ir^ii. For 
our simulations we chose a disc of 32 x 32 x 1 cells, with 
a grain size A of lOnm. 
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FIG. 4: Spontaneous equilibrium magnetization and parallel 
and transverse susceptibility vs. temperature. 



COMPARISON BETWEEN MULTI-MACRO 
SPIN MODEL AND EXPERIMENT 
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FIG. 5: Graph showing the electron temperature (Te) and the 
lattice temperature (Ti) simulated by the two-temperature 
model. Scaled reflectivity and magnetization data are also 
included for comparison. 



In order to make a comparison with the experimental 
data it is necessary to have an approximation to the tem- 
poral variation of the temperature pulse caused by the 
heating. Because of the complex behavior of the reflec- 
tivity as discussed by Djordjevic et al.— it is not feasible 
to use this property to determine the temperature profile 
in this work. Instead we make the simplifying assump- 
tion that, at low laser powers, the magnetization closely 
follows the electron temperature during the heating and 
initial recovery phase; an assumption essentially borne 
out by calculations in Refi^ using an atomistic model, 
where fast demagnetization and recovery were found for 
low peak electron temperatures. Following Ref,- , we as- 
sume that the photon energy is transfered to electrons 
and that the magnetization is directly coupled to the elec- 
tron temperature Te within a two-temperature model^S, 
expressed as the following coupled differential equations 
for the electron and phonon temperatures. To, T\ respec- 
tively: 



di; 

dTi 



-Gei(Te-Ti)+P(t) 



Ci--|^ = Gci{Tc — T\), 



(8) 



where Co, C\ are electron and lattice specific heat con- 
stant. Gal is a coupling constant and P{t) is the laser flu- 
ence. Eqs.[8]can easily be solved numerically to generate 
the time variation of Tc,T\. We determine the parame- 
ters of the two-temperature model by fitting to the form 
of the initial magnetization decay and recovery, assuming 
a Gaussian laser profile. 

Fig. [5] shows a comparison of the two-temperature 
model with the data from the reflectivity and low laser 
power magnetization measurements in order to estimate 
physically reasonable parameters for the model. The 
time-scale of the peak in electron temperature To matches 
the demagnetization peak in the magnetization data and 
the initial peak in the reflectivity data. The lattice tem- 
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FIG. 6: Calculated time dependence of the z-component 
of the magnetization for different applied magnetic fields. 

= 76 with K = 3.2 ■ 10^ J/m^ and V = (10 nm)^ at 
room temperature, iji = 1240K, A = 0.1) 



perature T\ reaches its peak value on the same time-scale 
as the second rise in the reflectivity data. For use in the 
computational model an interpolation function was fitted 
to the To predicted above. Different laser fluences were 
simulated by scaling the results to the peak electron tem- 
perature (TP), which becomes a main parameter in the 
comparison with experiment. 

We first describe calculations with the assumption of a 
spatially uniform temperature profile. The model param- 
eters used correspond to a material with a To of 660K, 
Ms — 1.75 xlO'^ A/m, an out-of-plane anisotropy K — 
3.2 X 10^ J/m'^ and a damping constant A of 0.1. The ma- 
terial is broken up into 32 x 32 cells of size 10 nm which 
are exchange decoupled in order to model the granular 
structure of a recording medium. The intention is to 
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FIG. 7: Calculated time dependence of the z-component 
of the magnetization for two different peak temperatures. 
^ = 76 with K = 3.2 • lO'^ J/m^ and V = (10 nm)^ at 
room temperature. (A = 0.1, Bz = — 0.52T) 
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FIG. 8: Calculated time dependence of the ^-component of 
the magnetization for different values of A. -^^^ = 76 with K 

= 3.2 • 10^ J/m^ and V = (10 nm)^ at room temperature {Ti 
= 1240 K, = -0.52T) 



outline the effects of the major parameters in the model, 
namely the applied field, peak electron temperature and 
the coupling parameter. Fig. [6] shows the temporal re- 
sponse of the magnetization to a laser pulse giving rise to 
a peak electron temperature of 12A0K assuming a cou- 
pling parameter of A = 0.1. It can be seen that the 
simulation gives a reasonable qualitative description of 
the time evolution of the magnetization following a laser 
pulse. Specifically we note the initial fast demagneti- 
zation and recovery. In the case of a positive field the 
magnetization recovers to the equilibrium value in the 
positive sense. In a reversing field the initial recovery is 
followed by a slow reversal of the magnetization toward 
the field direction which, in the model, can be unam- 
biguously attributed to thermally activated transitions 
over the particle energy barriers, supporting the earlier 
interpretation of the experimental data. 

Fig. [7] explores the effect of increasing the peak elec- 
tron temperature. As might be expected, the effect of 
this parameter is to lead to a more complete demagneti- 
zation during the laser pulse. After the high temperature 
pulse the system remains demagnetized due to the high 
temperature remaining after the pulse. 

Finally, we consider the effect of the coupling parame- 



ter A. Calculations for 3 different values of A are shown 
in Fig. [H Here, it can be seen that the effect of in- 
creasing A is to achieve a more rapid demagnetization. 
This is consistent with previous calculations^, where the 
effect is interpreted in terms of the more efficient trans- 
fer of energy into the spin system at large A leading to 
a more complete demagnetization for a given laser pulse 
width. Clearly the increased demagnetization caused by 
the stronger coupling to the conduction electron system 
results in an increased heat assistance of the magnetiza- 
tion reversal. 

In addition to the material parameters, the peak elec- 
tron temperature and the coupling constant are seen to 
be important in the heat assisted reversal process. Al- 
though the model reproduces the essential physics of the 
reversal process, particularly the different behavior on 
the timescales of longitudinal relaxation (fast demagne- 
tization) and transverse relaxation (super-paramagnetic 
fluctuations) qualitative agreement only is obtained with 
experimental data. Specifically, the value of KV/kBT — 
76 (at room temperature) used in the simulations is 
rather small in comparison with the experimentally de- 
termined values. In order to obtain closed agreement 
with experiment it is necessary to include an important 
experimental factor; specifically the fact that the probe 
pulse area is comparable to that of the pump. This in- 
troduces a significant temperature gradient within the 
probe area, which must be taken into account in the cal- 
culations. 

Here, we model this effect using a Gaussian tempera- 
ture profile for the laser spot. The diameter of the spot 
used in the model calculations is much smaller than in 
the experiments {800 fim). This would be an unrealis- 
tic approximation if the magnetization reversal involved 
domain wall processes; however, since the grains are es- 
sentially decoupled in the experimental films this sim- 
plified model is able to gives some insight into the ef- 
fects of the temperature profile. Since both pump and 
probe beams have a Gaussian temperature variation, we 
take the temperature variation to be of the form T(r) oc 
exp(— r^/rpjjjjjp), with rpump the radius of the pump 
beam. The MOKE signal is assumed to be proportional 
to the product of the magnetization with a sensitivity 
function. The sensitivity function is proportional to the 
area of material generating the MOKE signal at a partic- 
ular radius and the light intensity at that radius, which 
has a Gaussian weighting, i.e. cx nizr e'}qy{—r'^ /rpj.^y^^)dr , 
with r-probc the radius of the probe. A numerical integra- 
tion over the probe area was carried out to determine the 
calculated MOKE signal from the film. 

The introduction of a temperature profile results in a 
qualitative agreement with experiment using parameters 
close to the measured values. Fig. [S] shows calculations 
of the time dependence of the magnetization following a 
laser pulse with a peak temperature of 1480K for differ- 
ent values of the applied magnetic field. The parameters 
used were Aff = 0.4 x 10^ A/m, giving a room tem- 
perature value close to the measured (VSM) values of 
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0.32 X 10^ A/m. The anisotropy constant at zero tem- 
perature was taken as 3.94x lO^J/m'^, which gives a value 
of KV/k-QT = 95 at room temperature, in good agree- 
ment with the experimentally determined value of 96. It 
is interesting to note that if a larger value of is used 
then system appears to lock into a domain-state during 
reversal, which suggests that even at elevated tempera- 
tures the inter-granular magneto-static interactions can 
play an important role. The exchange coupling can be 
added, but values as big as 10% do not change the results 
significantly. As noted by Kazantseva et al.- , changes in 
the value of the damping constant affect the amount of 
energy the magnetic system absorbs from the initial heat 
pulse and so how much of a demagnetization is achieved 
for a given peak electron temperature. In addition there 
is a slight broadening in the demagnetization peak. 

Fig. [TOl shows a more detailed analysis of the LLB sim- 
ulations, indicating that the magnetization dynamics are 
an ultra-fast demagnetization and recovery caused by the 
electron temperature peak, after which the elevated tem- 
perature of the lattice causes a gradual switching of the 
individual grains of material. Fig. [10] shows the temporal 
evolution of the radial magnetization defined such that 
ma{r)dr is the spatially averaged magnetization over the 
annulus r — > r -t- dr; here a is the z-component of the 
magnetization to^ (Fig. llOf a)). the total magnetization 
(Fig. [TOfb)), and the in-plane magnetization, defined as 
the spatial average of (m^ -t- to^)^/^ (Fig. [roFc)). 

We consider the behavior of two distinct regions; the 
central region for radius < 7 cell radii, where the temper- 
ature exceeds during the pulse leading to complete de- 
magnetization, and the outer region which doesn't exceed 
Tc. The variation of the total magnetization (Fig. [TUf b)) 
is consistent with previous calculations^. Specifically, the 
rate of recovery of the magnetization depends upon the 
magnetic state. Within the central region the material is 
completely disordered and the recovery of the magneti- 
zation is relatively slow due to the need for the magne- 
tization to recover from highly disordered states. In the 
outer region the magnetization retains some memory of 
the initial state, which results in a rapid recovery^. 
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FIG. 10: Radial magnetization in a 0.52T field as a function 
of time after the laser pulse, (a) shows the perpendicular 
component niz, (b) the magnitude m, and (c) the magnitude 
of the in-plane component rrii 



Of most importance in terms of heat assisted reversal 
is the behavior of the central region. Heat assistance of 
the reversal is demonstrated clearly in Fig. fTOT a). which 
shows reversal of the central (heated) region while the 
magnetization in the outer region is not switched. The 
nature of the reversal in the central region is further in- 
vestigated using the radially resolved in-plane component 
of the magnetization, which is shown in Fig. [TOFc). It is 
interesting to note that a large in-plane component devel- 
ops on the timescale of 10 30 ps. This results from the 
relatively random recovery of the direction of the mag- 
netization after cooling though Tc. This contributes to 
the magnetization reversal in two ways. Firstly, some 
of the grains take on a negative sense of the magnetiza- 
tion on recovery. Others will recover in a positive sense 
but at an angle greater than the energy maximum as 
the anisotropy increases; these grains are most likely to 
switch into the negative direction as the temperature de- 
creases. At longer timescale, and consequently lower tem- 
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peratures, the in-plane component reduces as the magne- 
tization of each grain begins to he preferentially along the 
easy anisotropy axis. However, the in-plane component 
does not completely disappear, probably reflecting the 
Boltzmann distribution of the magnetization direction 
within each grain. In addition to these mechanisms it 
is also likely that there will be thermally activated rever- 
sal over the energy barriers at the elevated temperatures. 
Along with the increase in spontaneous magnetization as 
the film cools this would contribute to the gradual in- 
crease in niz over timescales of 1 ns. 

VI. CONCLUSIONS 

We have presented laser pump-probe measurements 
which show a clear heat assistance from the laser pulse for 
switching the magnetization state. This is demonstrated 
by the ability to switch in an externally applied field with 
a magnitude lower than the intrinsic coercivity. The ex- 
periments show a rapid demagnetization and recovery fol- 
lowed by a slow evolution of the magnetization into the 
field direction. This is consistent with the existence of 
two characteristic relaxation times; the longitudinal and 
transverse relaxation times. The former is atomic-scale 
processes and is typically of the sub-picosecond order, 
whereas the transverse relaxation time reflects transitions 
over the energy barrier and can be orders of magnitude 
longer. In order to investigate the reversal mechanism we 
have developed a micromagnetic model based on the LLB 
equation which naturally includes both timescales. The 



LLB model calculations are in good quantitative agree- 
ment with the experimental data, as long as the the tem- 
perature gradient across the probe pulse is included. It 
appears that on the short time-scale (2 ps) there is a rapid 
demagnetization of due to an associated loss of Ms 
via the longitudinal relaxation. There is a partial recov- 
ery oiruz in the original direction as Ms starts to recover. 
However, the switching is assisted by the recovery of the 
magnetization of individual grains in random directions 
as the system cools through Tc . On the longer time-scale 
the reversal of in the applied field may also be assisted 
by thermally activated switching caused by the elevated 
lattice temperature. The elevated temperature has the 
effect of lowering the anisotropy energy barriers (due to 
the reduced values of the Magnetocrystalline anisotropy 
energy) and also provides the thermal energy to induce 
the transitions. The complex behaviour requires a model 
including both the longitudinal and transverse relaxation 
times, which is included here using the LLB equation. 
Our LLB based calculations encapsulate the physics of 
the heat-assisted reversal process and suggest the LLB 
equation as a physically realistic model for Heat Assisted 
Magnetic Recording. 
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